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. Abstract. For a positive integer n, let Mn be the set oi nxn complex matrices. Suppose m, n > 2 

are positive integers and k £ {1, . . . , mn — 1}. Denote by WkiX) the fc-numerical range of a matrix 
X £ Mmn- It is shown that a linear map <j) : Mmn ~^ Mmn satisfies 

; Wk{(p{A(X)B)) = Wk{A(X)B) 

1/^ ' for all A £ Mm and B £ M„ if and only if there is a unitary U £ Mmn such that one of the following 

holds. 

(i) For all A G A/„, B £ M„, 

(P(A»B) = Ui^iA»B))U*. 

(ii) mn = 2k and for all A £ Mm, B £ M„, 

(j>{A®B) = (tr {A ® B)/k)Imn - U(if{A ® B))U' , 

^ , where (1) (p is the identity map A0B^A®B or the transposition map A0 B ^ (A® BY , or 

(2) min{m, n} < 2 and <p has the form A(8)BM>A®B'orA(g)Bi-^A'®B. 

2010 Math. Subj. Class.: 15A69, 15A86, 15A60, 47A12. 
' i^'ey words: Hermitian matrix, linear preserver, A;-numerical range, tensor product of matrices. 
I> . 

1. Introduction and the main theorem 

CN . Let Mn {Hn) be the set of nxn complex (Hermitian) matrices. For k E {1,...,?7-}, the k- 

■ numerical range of A € Mn is defined as 

Wk{A) = {tr {X*AX)/k : X isnxk, X* X = h}. 

^ ' Since Wn{A) = {tr(j4)/n}, we always assume that k < n to avoid a trivial consideration. When 
A; = 1, we have the classical numerical range VFi(^) = VF(^), which has been studied extensively; 
see |10 1 I1 H [T2]. Denote hj A^B the tensor (Kronecker) product of matrices A € Mm and B € M^. 
The purpose of this paper is to characterize linear maps on M^n satisfying 

(1) Wk{A ®B) = Wk{(^{A ® B)) for ah A G M^, B e M„. 

The study is motivated by two areas of research. 

First, in the last few decades there has been considerable interest in studying linear preservers 
on matrix algebras as well as on more general rings and operator algebras. For example, Frobenius 
[8] showed that a linear operator (j) : Mn Mn satisfies 

det(</.(^)) = det(^) 

for all A € Mn if and only if there are U,V Mn with det{UV) = 1 such that <j) has the form 

(2) A^UAV or A^UA^V, 

1 



2 AJDA FOSNER, ZEJUN HUANG, CHI-KWONG LI, YIU-TUNG POON, AND NUNG-SING SZE 

where denotes the transpose of A. Clearly, a map of the above form is linear and leaves the 
determinant function invariant. But it is interesting that a linear map preserving the determinant 
function must be of this form. Furthermore, in [3] Dieudonne showed that a bijective linear operator 
(j) ■ Mn — )• Mn maps the set of singular matrices into itself if and only if there are invertible 
U,V £ Mn such that 4> has the standard form ([2|) . One may see [16] and references therein for 
results on linear preserver problems. For more new directions and active research on preserver 
problems motivated by theory and applications, one may see, for example, [Tj [20l [25]. 

In connection to the linear preservers of the /c-numerical range, Pellegrini [22] proved that every 
linear map (p '■ — > Mn preserving the numerical range is of the form 

(3) A ^ UAU or A^ UA^U* 

for some unitary U € M„. Three years later. Pierce and Watkins [23] extended this result to the 
fc-numerical ranges as long as n 7^ 2k. In [1^ (see also [21]) it was shown that for n = 2k, a linear 
map (f) : Mn — > Mn preserves the /c-numerical range if and only if there exists a unitary U G M„ 
such that (j) has the form ([3|), or 

(4) A ^ (tr {A)/k)In - UAU* or A ^ (tr {A)/k)In - UA*U*. 

One may see [21 for more information about the results on linear maps on M„ which preserve 
the fc-numerical range. 

Another motivation of our study comes from quantum information science. In quantum physics 
(e.g., see [9]), quantum states are represented by density matrices D in Hn, i.e., positive semidefinite 
matrices with trace 1. If D has rank one, i.e., D = xx* for some x € C" with x*x = 1, then D is a 
pure state. Otherwise, D is a mixed state, which can be written as a convex combination of pure 
states. In a quantum system, every observable corresponds to a Hermitian matrix A. Then 

W{A) = {tr (Axx*) : X G C", x*x = 1} 

can be viewed as the set of all possible mean measurements of quantum states. If ^ = ^1 + iA2 € 
Mn, where Ai,A2 S Hn, then every point in VK(^) and be viewed as the set of all joint measurements 
x*Ax = x*Aix + ix* A2X of the quantum state x with respect to the two observables associated 
with Ai and A2, and thus VF(74) is the set of all possible joint measurements. By the convexity of 
W{A) (e.g., see [El dH HI] ) , 

W{A) = {tr {AD) : D e Hn is a density matrix}. 

Thus, 14^(^) can also be viewed as the set of all possible joint measurements of two observables 
on mixed states. Now, by the convexity of Wk{A), and the fact that the convex hull of the set 
{P/k : = P = P* ,ti P = k} equals the set Sk of density matrices D satisfying In/k — D is 
positive semidefinitie, we have 

Wk{A) = {tr{AP)/k : P € Hn,P^ = P = P* ,tT P = k} = {ti {AD) : D € Sk}. 

So, VFfc(74) can be viewed as the set of joint measurements of the states in Sk with respect to the 
observables associated with A. Suppose A £ Hm and B G Hn correspond to observables of two 
quantum systems with quantum states Di G Mm and D2 ® Mn. Then the tensor (Kronecker) 
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product A (S) B correspond to an observable on the joint (bipartite) system with quantum states 
D G Mmn expressed as the convex combination of uncorrelated quantum states Di (8) D2, where 
Di G Mm, and D2 G M„ are quantum states (density matrices). In this connection, we are interested 
in studying hnear maps (j) on Mmn satisfying ([1]). 

In fact, this hne of study has been carried out in [5l [H [71 [131 [2l]- Suppose f{X) is a hnear 
function on the matrix X G Mmn- It is shown in [H] that the hnear maps (p on Hmn satisfying 

(5) f{ct){A(^B)) = f{A(^B) 

for ah A G Hm and B G Hn when f{X) is the spectrum or the spectral radius oi X . In [6l [7], the 
authors characterized the linear maps 4> on Mmn satisfying ([5|) for all A G Mm and B G Mn when 
f{X) is the Ky Fan norms, Schatten norms or the numerical radius of X. 

The following is our main result. 

Theorem 1.1. Let k G {1, . . . , mn — 1}. ^ linear map (j) : Mmn Mmn satisfies 

(6) Wk{HA(^B)) = Wk{A(^B) 

for all A G Mm and B G M„ if and only if there is a unitary U G Mmn such that one of the 
following holds. 

(i) For all A G Mm, B G M„, 

(7) (l){A(^B) = U{f{A(^B))U*. 

(ii) mn = 2k and for all A G Mm, B G Mn, 

(8) (P{A ^B) = (tr (.4 B)/k)Imn - U{<f{A ® B))U* , 

where (1) ip is the identity map A B A ^ B or the transposition map A0 B >—?■ {A ^ B)* , or 
(2) min{m, n} < 2 and ip has the form A(^B 1-^ A(^B^ or A B 1-^ A^ B . 

The proof of the theorem will be given in the next section. We will use the following properties 
of the /c-numerical range; for example, see [11 1 114^ [T9l I23j. 

Proposition 1.2. Let A G M„ and A; G {1, . . . , n — 1}. 

• For any a, /3 G C, Wk{A) = {a} if and only if A = aLn, and Wk{aLn + (3A) = a + /3WkiA). 

• For any unitary U G Mn, WkiUAU*) = Wk{A). 

• For any s x n matrix V with s > k and VV* = Lg, we have Wk(VAV*) C Wk{A). 

• VFfc(^) QM if and only if A is Hermitian. 

• If A £ Hn has eigenvalues «!>•••> an, then 

WkiA) = [{on-k+i + ■■■ + an)/k, (tti + • • • + ak)/k]. 

• Wk{^) = Re {Wk{A)) = {Re (z) : z G Wk{A)}. 
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2. Proof of Theorem 11.11 

In the following, denote by Eij G M„, 1 < i,j < n, the matrix whose («, j)-entry is equal to one 
and all the others are equal to zero. The ^-numerical radius of A S M„ is defined as 

Wk{A) = max{\x\ : x G Wk{A)}. 

Two matrices A,B M„ are called orthogonal if AB* = A*B = 0. We write AJ-B to indicate that 
A and B are orthogonal. It is shown in [18] that A^B if and only if there are unitary matrices 
U,V e Mn such that UAV = Diag (ai, . . . , and UBV = Diag . . . , /3„) with ai,/3i > 
and aj/3j = for i = l,...,n. The matrices ^i,...,^^ are said to be pairwise orthogonal if 
A*Aj = AiAj = for any distinct i,j G {l,...,s}. In this case, there are unitary matrices 
U,V Mn such that UAiV = Di for i = 1, . . . , s with each Di being nonnegative diagonal matrix 
and DiDj = for any distinct i,j G {1, . . . , s}. We will need the following lemmas in the proof. 
The first lemma was proved in |15j . 

Lemma 2.1. |15| Lemma 4.1] Let k G {l,...,n} and suppose A G have diagonal entries 
ai,...,a„ and eigenvalues Ai > • • • > A„ respectively. Then Ylj=i'^j ~ Sj=i -^j ^/ ^''^^ ^^^^2/ ^/ 
A = Ai (B A2 where A\ G B]^ has eigenvalues Ai, . . . , A^. 

Lemma 2.2. Let k G {l,...,n} and A,B£ Hn he positive semidefinite matrices. Suppose 
tT{A)/k = max{x : x e Wk{A - B)} . Then A±B. 

Proof. Suppose U G M„ is a unitary matrix such that U{A — B)U* = Diag (Ai, . . . , A„) with 
Ai > ••• > A„. Denote the diagonal entries of UAU* and UBU* by ai,...,a.„ and 61,..., 6„, 
respectively. Then a,,, 6j are nonnegative for i = 1, . . . ,n, since ^ > and B > 0. Now 

k k 

Y^{ai-bi) = ^Xi = tr (A) 

i=l i=l 

leads to 

k 

^ aj = tr (A), Ofc+i = • • • = a„ = and 61 = • • • = 6^ = 0. 

1=1 

Using the fact A>0 and B >0 again, UAU* and UBU* must have the form 

UAU* = Ai(B On-k, UBU* = Ofc e -Bi 
which means AJ-B. □ 

Denote by Xi{X) > • • • > A„(X) the eigenvalues of a Hermitian matrix X G M„. 

Lemma 2.3. Lei A; G {1, . . . ,mn — 1} anrf (p : Mmn M^n be a linear map satisfying The 
following conditions hold. 

(a) (j){Hmn) ^ Hmn- 

(b) (t>{Lmn) = Imn- 

(c) (/> is trace-preserving. Furthermore, if mn 7^ 2k, then (pi^Ea (g) is positive semidefinite 
for all 1 < i < m and 1 < j < n. 
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Proof, (a) Let A £ Hm and B € H,,. Then Wk{(l){A (g, B)) = Wk{A (g)B) CR. By Proposition 
fOl 0(A ®B)CR. Since every C S -f^mn is a linear combination of matrices of the form A (g) B 
with A G and B G we see that (j) maps Hmn to Hmn- 

(b) W^fc(</.(/™„)) = Wu{In.n) = {1}. Thus, = Ir,^n■ 

(c) Let ai > • • • > be eigenvalues of ® Ejj) = Aij. Since Wk{Aij) = WkiE^ (g) E'jj) = 
[0, l/k], we have ai + - ■ ■ + ak = 1 and am„_fc+i + • • ■ + amn = 0. If mn = 2k, then tr Aij = 1 + = 1. 

If mn > 2k, then a^+i, . . . , a^n-fc ^ oimn-k+i ^ and, thus, tr(^jj) > 1. Moreover, if 
otmn-k > 0, then tr (Aij) > 1. On the other hand, we have 

mn = tr {Imn) = tr {(pilmn)) = tr ^'/'(^^ ® -^'ji) j = tr Ajj^ > mn. 

This yields that Aij is a positive semidefinite matrix with trace one. 

Similarly, if mn < 2k, then amn-k+i + ■ ■ ■ + Ofc ^0 and, thus, tr {Aij) < 1. Therefore, 

mn = tr {Ir,^n) = tr {(j){Imn)) = tr ^(Yl -^ii)^ = - 

which yields that Ay is a positive semidefinite matrix with trace one. 

We can apply the same argument to show that for any orthonormal bases {xi, . . . ,x„i} C 
and {ui, . . . , Un} C C^, tr {(j){xiX* (8) UjUj)) = 1- Thus, (p is trace preserving for all Hermitian A^ B 
and, hence, for all matrices in Mmn- D 

Lemma 2.4. Let k € {2,...,N -2} and X,Y e Hn with Wk{X) = WkiY) = [0,1/k] and 
WkiX + Y) = [0,2/k]. If 

X -Y = Diag (1 - (A; - l)o, a, . . . , a, -1 - (fc - l)a) 

with a € [— l//c, 1/k], then 

X = Diag (1 - (A; - l)d, d,...,d,-{k- l)d) and Y = Diag {-{k - l)d, d, . . . ,d,l - {k - l)d) 

with d G {0, 1/A;} so that X - Y = Diag (1,0,... ,0,-1). 

Proof. Suppose X has diagonal entries xi, . . . , xj^f and Y has diagonal entries yi, . . . , un- Then 
for any 1 = ii < . . . < < iV — 1, we have 

k k k k 

^Xi^<l, ^yit>^, ^Xi^-^yi^ = l, 

i=l t=l t=l t=l 

which imply Xj^ + • • •+Xi^ = 1 equal to the sum of the k largest eigenvalues of X and + • • •+yi^. = 
equal to the sum of the k smallest eigenvalues of Y. Thus, applying Lemma |2. II Lemma 4.1]), 
xi is the largest eigenvalue of X, X2 = • • • = x^v-i are the second largest eigenvalue of X, yi is the 
smallest eigenvalue of Y, and y2 = ■ ■ ■ = Vn-i are the second smallest eigenvalue of Y. Moreover, 
X and Y have the form 

X = Diag(l -{k- l)x,x,.. .,x,-{k - l)x), Y = Diag (-(A; - l)y,y, ...,y,l-{k- l)y) 
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with x,y £ [0, Hence, X+Y = Diag {l — {k—l)x — {k — l)y, x+y, . . . , x+y, l — {k—l)x — {k—l)y) 
satisfies Wk{X + Y) = [0, 2/A;]. So, either (a) k{x + y) = 2, which imphes that a; = y = 1/A;, or (b) 
k{x + y) = 0, which imphes that x = y = 0. □ 

Lemma 2.5. Let 2 < k < mn/2 he an integer and (j) : Mmn be a linear map satisfying 

Then for any orthonormal bases {xi, . . . , Xm} Q C™ and {yi, . . . , yn} ^ C", either 

(1) there is a unitary U £ Mmn such that U*(j){xiX*®yjy*)U = XiX* ^y-jy* for all i = 1, . . . ,m, 
and j = I, . . . ,n, or 

(2) mn = 2k and there is a unitary U G Mmn such that U* (j){xiX*®yjy*)U = Imn/k—XiX*0yjy* 
for i = 1, . . . ,m, and j = 1, . . . ,n. 

The proof of this lemma is rather technical, we will present it in the last part of this paper. 
Denote by ^{X) the spectrum of X € Mn- The following example is useful in our proof. 



Example 2.6. Suppose m,n> 3. Let A = X ® Om-3 and B = X ® On-3 with X 
Then A® B is unitarily similar to 
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and A0 B^ is unitarily similar to 
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Consequently, 



a{{A(^B + iA^B)*)/2) = {-741/2, -3/2, -3/2, 0, 0, 3/2, 3/2, V41/2}, 
a{{A(^B* + {A(^B^)*)/2) = {-9/2, -^972, -1/2, 0, . . . , 0, 1/2, ^9/2, 9/2}. 

Applying Proposition [Ol one see that Re {WkiA(g)B)) / ReWk{{A®B^)), and hence Wk{A(^B) 7^ 
Wk{A (g) i?*) for any k e {1, . . . , mn - 1}. 



Now we are ready to present the proof for Theorem II. li 
Proof for Theorem II. li 

Note that the 1-numerical range is just the classical numerical range. The case k = 1 has been 
obtained in [71 Theorem 2.1]. Since (n — k)Wn-kiA) = tr (A) — kWk{A), by Lemma [231 we have 

Wki^A)) = Wk{A) ^ Wn^kiH^)) = Wn-k{A) . 

Therefore, we can focus our proof on 2 < A; < mn/2, with mn > 4. 

Note that Wk{X) = Wk{^{X)), Wk{X) = Wk{X^), Wk{X) = Wk{U*XU) for any unitary U 
and X in Mmn- Furthermore, if ^4 € M2, then A = UaA^Ua for some unitary Ua depending 
on A so that for any B G M„, the matrix A 1^ B is unitarily similar to A* (8) -B. Thus, Wk{A (8) 
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B) = Wk{A* ^ B) = WkiiA^ ® BY) = WkiA O S*). Similarly, if ^ G M„ and B € Mg, then 
Wk{A (g) i?) = Wk{A (g) i?*) = WkiA^ ® i?). Combining the above, we get the sufficiency. 

For the converse, suppose WkiA B) = Wfc(</)(A (g) B)) for all {A,B) £ Mm x M„. Suppose 
mn ^ 2k. Then by Lemma |2.5|, (1) always holds. So, for any Hermitian A E Mm and B € M„ 
with spectral decomposition A = ^ajXjX* and B = X^&j^/jyj, we see that (j){A (g) B) is unitarily 
similar to ACS)B. So, y4(gi? and (j)(A(^B) always have the same eigenvalues. Thus, by [H Theorem 
3.2], there is a unitary V such that (p has the form A ® B V*ip{A (g B)V for any Hermitian 
A G Mm and i3 G M„, where (/3 is one of the following maps: A(giB\-^A®B,A®B\-^A(i^ B^, 
A®B^A^^BoxA®B^A^®B^. By linearity, the map (j) can only have one of these forms 
on Mmn- However, if m, n > 3, we see that ip cannot be of the second and third form by Example 
12.61 So, if can only be of the first and fourth form. The desired conclusion holds. 

Now, suppose mn = 2k. We claim that either (1) in Lemma [2 . 5 1 holds for any choice of orthonor- 
mal bases {xi, . . . , Xm} Q C™ and {yi, . . . , y„} C C", or (2) in Lemma 12.51 holds for any choice of 
orthonormal bases {xi, . . . , Xm} ^ C™ and {yi, . . . , y„} C C". To see this, note that the set 

S = {(x, y) : X G C™, y G C", x*x = l = y*y} 

is path connected because the unit spheres in C™ and C" are path connected. Consider the 
continuous map from 5 to reals defined by f{x,y) i— > | det(0(xx* (g yy*))\- If (1) holds for a pair 
of orthonormal bases containing x and y, then f{x,y) = 0; if (2) holds for a pair of orthonormal 
bases containing X and y, then /(x, y) = | det(Imn/^ — ^ii (^i-Eii)! = (1/A;)'""~^(1 — l/Zc). It follows 
that either /(x, y) = for all (x, y) G S so that (1) always holds, or /(x, y) = — 1/k) for 

all (x, y) G S so that (2) always holds. 

If (1) holds for all orthonormal bases {xi, . . . ,Xm} ^ C'" and {yi, . . . Q C", then by the 
argument as in the case of mn ^ 2k, we see that (p has the desired form. If (2) holds for all 
orthonormal bases {xi,...,Xm} Q C™" and {yi,---,yn} Q C", then compose (j) with the map 
X 1-^ (tr X)I/k — X so that the resulting map satisfying (1). The result follows. □ 

Finally, we give the the proof of Lemma 12.51 
Proof for Lemma 12.51 

By Lemma 12.31 (a), (p maps Hermitian matrices to Hermitian matrices. We may focus on the 
case that XjX* = En for i = 1, . . . ,m and yjy^ = Ejj for j = 1, . . . ,n. Otherwise, replace <j) by the 
map (f){A'S) B) = (piViAV^ (g V2-BV^2*) s° ^^at Vi G Mm and V2 G Mn are unitary matrices satisfying 
ViXiX*V{ = Ejj for i = 1, . . . , m, and V2yjy*V2 = Ejj for j = 1, . . . , n. 

We divide the proof into three cases, namely 

(a) mn 7^ 2k, (b) mn = 2k, m < 3 and n < 3 and (c) mn = 2k, m > 4 or n > 4. 
2.1. The case mn ^ 2k. 

Claim 1 . There exists a unitary U G Mmn such that 

cP{Eii g) Ejj) = U{Eii g) Ejj)U* 

for all 1 < i < 77T, and 1 < j < n. 
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It suffices to prove that 

(9) (/){Eii ® Ejj)±^{Err E,s) 

for all pairs {i, j) ^ (r, s) with 1 < i,r < m and 1 < j, s < n. 
First, suppose that i = r or j = s. Considering 

Wk{(p{E,, (g) Ejj) - <j){E„ (8) E,s)) = Wk{Eii <S) Ejj - E^ E.^) = [-l/k, 1/k], 

applying Lemma 12.31 and Lemma 12.21 we conclude that ([9]) holds. 

Now, let i ^ r and j ^ s. We may assume that 2 < k < mn/2 < mn — 2, we consider 

Wk{(l){EiMEjj+Ess))-HErMEjj+Ess))) = Wk{EiMEjj+Ess)-ErrMEjj+E,s)) = [-2/k,2/k]. 

Applying Lemma 12.21 again, it follows that <j){Eii (8) {Ejj + Ess))-^(t>{Err ® (Ejj + -Ess))- Hence, we 
have dH). 

2.2. The case mn = 2k, m < 3 and n < 3. 

Since mn = 2k is an even integer, without loss of generality we may assume that n is even. So 
it suffices to consider the cases when n = 2 and m G {2, 3}. 

Claim 2 . Let m G {2, 3} and Ai = (j){Eii ® [En - E22)) € Mm M2 for f = 1, . . . , m. Then there 
is a unitary U G such that (j){Ai) = U {En (g) {Eu — E22))U* for i = 1, . . . , m. 

We need only show that Ai,...,Am are mutually orthogonal and each Ai has eigenvalues 

1, -1,0, ... ,0. Note that Wk{Ai) = [-l/k, l/k] for i = 1, . . . , m. So Ai(^i) H h XkiA) = 1 and 

Afc+i(Ai) + • • • + X2miAi) = -1. Since WkiAi + A2) = [-2/k, 2/k], we see that 

k k 

X,{Ai + ^2) = J2(Xj{A^) + A,(^2)) = 2 

and 

2m 2m 

X^{A, + A2)= Y (Ai(Ai) + A,(^2)) = -2. 
j=k+i j=k+i 

So, by a unitary similarity and applying Lemma l2.ll we may assume that Ai = Bi (B Ci and A2 = 
B2®C2 so that Bi has eigenvalues \i{Ai), . . . , \k{Ai) and Ci has eigenvalues Xk+i{Ai), . . . , X2m{Ai), 
i = 1,2. As WkiAi - A2) = [-2/k, 2/k], we see that {Bi - B2) © (Ci - C2) has eigenvalues 

71 > ■ • ■ > 72m such that 71 H h 7fc = 2 and 7^+1 H h 72m = -2. Clearly, 71, . . . , 7^ cannot 

all come from Bi — B2, else, 71 + • • • + 7fc = tr {Bi — B2) = 0. Similarly, 71, ■ . ■ , 7fc cannot all come 
from Ci — C2. Now we distinguish two cases. 

Case 1. m = 2. In this case we see that an eigenvalue of Bi — B2 and an eigenvalue of Ci — C2 
sum up to 71 + 72 = 2. Since \i{Bi — B2) < Ai(yli) — A2(^2) and 



Al(Ci-C2) < A3(Ai)-A4(A2), 
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we have 

2 = 7l+72<Ai(Ai)+A3(Al)-A2(A2)-A4(A2) 

< Xi{Ai) + X2iAi) - X3{A2) - Xa{A2) = 2. 

It follows that 

A2(^i) = A3(^i) and X2{A2) = X3{A2). 
Without loss of generality, assume that Ai is unitarily similar to a matrix of the form 

Diag(l -a,a,a, -1 -a) with a G [-1/2, 1/2]. 

By Lemma |2.4|, we conclude that a = 0. 

Similarly, we can show that A2 has eigenvalues 1,-1,0,0. It is then easy to show that Ai,A2 
are orthogonal. 

Case 2. m = 3. We have two subcases. 

Subcase 2.1. An eigenvalue of Bi — B2 and two eigenvalues of Ci — C2 sum up to 71+72 + 73 = 2. 
Since Xi{Bi - B2) < Xi{Ai) - X2,{A2) and 

Xl{Ci-C2) + X2{Ci-C2) < Ai(Ci) + A2(Ci)-A2(C2)-A3(C2) 

= A4(^l) + A5(Ai)-A5(A2)-A6(A2), 

we have 

2 = 7l+72 + 73 <Ai(^l) + A4(^l) + A5(^l)-A3(^2)-A5(^2)-A6(A2) 
< Ai(Ai) + X2{Ai) + XsiAi) - A4(A2) - A5(A2) - A6(^2) = 2. 

It follows that 

A2(^i) = --- = A5(Ai) and X3{A2) = Xi{A2). 

Subcase 2.2. An eigenvalue of Ci — C2 and two eigenvalues of Bi — B2 sum up to 71 + 72 + 73 = 2. 
Then 

A3(^i) = A4(^i) and A2(^2) = • • • = A5(^2). 

Without loss of generality, assume that Ai is unitarily similar to a matrix of the form 

Diag (1 — 2a, a, a, a, a, — 1 — 2a) with a G [—1/3, 1/3]. 

By Lemma 12.41 we conclude that a = 0. 

Now, applying the arguments to A2 and A^, we conclude that one of the matrices A2 and A3, say, 
A2, has eigenvalues 1, —1, 0, 0, 0, 0. Then it follows from Lemma [2.2l that Ai and A2 are orthogonal. 
So, we may assume that Ai = Diag (1,0, 0, —1,0, 0) and A2 = Diag (0, 1, 0, 0, —1, 0). Note that 

WsiAi +A2- A3) = W^iAi +A2 + A3) = [-1, !]■ 

We see that ^3 = i?3©C3 such that i?3 has eigenvalues 1— 2c, c, c and C3 has eigenvalues c, c, — 1 — 2c. 
Now, applying Lemma [2^ on A3, we conclude that A3 also has eigenvalues 1, —1, 0, 0, 0, 0. It follows 
from Lemma 12.21 that Ai, A2, A3 are mutually orthogonal. Thus, we obtain the claim. 
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Using the notation as in Claim 2, we see that there is a unitary U € such that Ai = 

U{Eii (g) {Ell - E22))U*. By Lemma EH for each i = 1, . . . ,m, 

(10) <j){Eii ® Ejj) = U{Eii (g) Ejj)U\ j = 1, 2, 
or 

(11) (l,{Eu Ejj) = hm/k - Pi{Eu ® Ejj)Pl, j = 1, 2, 

for a suitable permutation matrix Pi € M2m- Because Yl,ij4>{Eii ® Ejj) = 12m.-, either ([TO]) holds 
for all i = l,...,m, or (llip holds for all z = 1, . . . , m. In the latter case, the map 4'{A) = 
tr {A)/kImn — 4>{-^) must satisfy (fTOl) . This shows that Pi = • • • = Pm- 

2.3. The case mn = 2k, m > 4 or n > 4. 

Without loss of generality, we assume m > 4. If Aij = (p{Eii (g) Ejj) is positive semidefinite for 
any 1 < i < m and 1 < J < then applying the same arguments as in the previous case, we 
conclude that (p satisfies (1). Now suppose there exists some iq and jo such that Ai^j^^ has negative 
eigenvalues. Without loss of generality, we assume io = jo = 1- 

Claim 3 . There exists some i G {1, . . . , m} such that (^{Eii ® In) has a negative eigenvalue. 

For 1 < i < m, we denote the eigenvalues of (^[Eu ® In) by ai{i) > a2{i) > • • • > amn{i)- Since 
Wk{4>{Eii eg In)) = [0, n/k]^ we have Yl^j=i %(*) = ™d Yl^=k+i ^ji''') = for 1 < i < m. Suppose 
^{Eii^In) > for all 1 < i < m. Since Wk{(t){Eii+E22)®In) = [0,2?i//c], without loss of generality 
we can assume 4>{{Eii + E22) ig /«) = Diag (ri, . . . , r^n) with ri > r2 > ■ ■ ■ > rj, > r^+i = ■■■ = 
Tmn = 0. Let (t){Eii (g /„) = {xij) and (/)(£'22 <g In) = (Uij)- Then 

k k k k 

'^Xii = ^yii = n = ^ai(l) = ^0^(2). 

i=l i=l 1=1 1=1 

By Lemma [27T1 (/"(-E^ii (g /„) = X © 0^ and 4^{E22 ® In) = Y ®Qk with tr X = tr T = n. Moreover, 
W}^{(j){{Eii — E22) (g In)) = Thus, applying Lemma 12.21 we have X1.Y . It follows that 

X is singular and afc(l) = afc+i(l) = • • • = am„(l) = 0. Suppose V G -Mmn is a unitary matrix such 
that 

n 

Diag (ai(l), . . . , afc_i(l), 0, . . . , 0) = ^.^(^11 ® In)V* = y,^(£;ii ® Ejj)V*. 
Denote the diagonal entries of V(j){Eii Ejj)V* by di(j), . . . , dmnij)- Then 

fe— 1 n fe— 1 

n = ^ a,(l) + a,(l) = ^(^ d,{j) + d,(j)) 

i=l j=l i=l 

for every s € {k, . . . , mn} and VFfc((^(ii^ii -Ejj)) = [0, 1/A;] ensures that J2i=i di{j) + ds{j) = 1 for 
every s € {/c, . . . , mn}. Applying Lemma |2. II again, we see that for every j € {1, . . . , n}, we have 

V()){Eii (g Ejj)V* = Rj © tjik+i, 
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where each eigenvalue of Rj is larger than or equal to tj. Further, G Wk{4>{Eii Ejj)) implies 
tj = for 1 < j < n, which contradicts with the assumption that (j){Eii -En) is not positive 
semidefinite. 

Claim 4- Suppose there exists i G {1, . . . ,m} such that the eigenvalues of (j){Eii (X" In) are ai{i) > 

O-mnii) with (Xmn 

(i) < 0. Then 

(12) ai{{} = a2{i) = ■ ■ ■ = ak+iii). 
Moreover, (piEjj (8) In) has a negative eigenvalue for every j G {1, . . . , m}. 

Without loss of generality, we assume i = 1 and 

(p{En (g) /„) = Diag (ai(l), • • • ,amn(l)) 

with 

(13) ai{l)>a2{l)>--->amn{l), 

where ai(l), . . . ,afc_|_i(l) are not identical. Let us denote the diagonal entries of (piEjj (8) In) by 
hi{j), . . . , hmn{j) and the diagonal entries of U(p{Ejj (g) In)U* by hi{U,j), . . . , hmniU,j). Note that 
afc(l) and ak+i{l) must be equal. Otherwise, by the fact that Wk{(p{{Eii + Ejj) ^ In)) = [0,2n/k], 
we have Ylt=i "■r{'^) = Ylr=i ^r{j) = n for j = 2, . . . , m. But then, if Z = Y^f^i (p{Ejj (g) /„), 
the leading k x k submatrix of Z will have trace mn, which contradicts with Wk{4>{Imn)) = {I}- 
Suppose ai(l) > afc(l)i i-e., there are integers s,t G {1, . . . , k — 1} such that 

(14) ai(l) > • • • > as{l) > as+i{l) = ■■■ = afc+t(l) > afc+t+i(l) > • • • > a„„(l). 
We are going to show that 

(15) hiij) = ■■■ = hs{j) = hk+t+iU) = ■■■ = hmnU) for j = 2, . . . , in. 
Let 7 = m/2 when m is even and 7 = (m + l)/2 when m is odd. Denote by 

G = (j)i{2{Eu + ■■■ + ^7-1,7-1) + ^7,7) <S) In), 

Gi = <^((^ii + ••• + £;^_i,^_i) ®I„), and G2 = (^((Sn + • • • + ^^,^) I„). 
Then we have 

(16) WkiGi) = [0, (7 - l)n/A:], Wk{G) = Wk{Gi) + Wk{G2), 

where the fc-numerical radius Wk{G), Wk{Gi), and Wk{G2) are the right end point of Wk{G), Wk{Gi), 
and Wk{G2), respectively. Let ^7 be a unitary such that the sum of the first k diagonal entries of 
UGU* equals to kwk{G). Then the sum of the first k diagonal entries of UGiU* equals to kwk{Gi) 
for i = 1,2. We assert that the following conditions hold. 

(a) Ylp=i ^pi^il) — when m is even and '}2^=i ^p(^)7) = when m is odd. 

(b) U(i){Ejj ® In)U* = Bji © Bj2 with Bji € and tr (Bji) = ?i for j = 1, . . . , 7 - 1. 

(c) U(j){Ejj (g) In)U* = Bji © Bj2 with Bji G Mk and tr {Bj2) = n for j = 7 + 1, . . . , m. 
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Since 

(7 - l)n = kwkiGi) = Y.U Eti hr{U,j) < EU n = {l-l)n, and 

— = kWk{G2) = Ei=i Er=l hr{U,3) < —. 

It follows that X]p=i ^p(f^5 j) = ^ i = • • • ; 7 ~ 1 {^) holds. Applying Lemma |2.H we have 
the condition (b). 

For any j G {7 + 1, . . . , m}, since 

WkWEjj In)) = Wk{Ejj /„) = [0, n/k], 

the sum of any k diagonal entries of U(f){Ejj (8) In)U* lies in [0, n]. Now, the right end point of the 
set 

Wk{<p{{Eii + • • • + + Ejj) ® /„)) 
is 1, and the sum of the first k diagonal entries of U4>{(Eii + • • • + -£^7,7) (8) In)U* is k. We see that 
Ep=i hp{U,j) = and E^fc+i hp{U,j) = n. Hence, we get (c). 
Suppose Ui and U2 are unitary matrices such that 

C/iBiif/i* = Diag(ai(l),...,afc(l)) and [/a^iaf/l = Diag (afc+i(l), . . . , a^„(l)). 

Replace U with (C/i©C/2)f^- Then the new matrix U also satisfies (a), (b) and (c). Moreover, U is of 
diagonal block form i/s © C/4 f/5 with G M^, f/s G Mf^^t. Since any unitary C/3 and C/5 will yield 
the same summation of the first k diagonal entries of UGU* , we can assume f/ = /g C/4 Ik-t- 
Thus, for j = 7 + 1, . . . , m, it follows from (c) that 

hp{j) = hp{U,j) <hq{U,j) = hq{i) for aU p G {1, . . . ,s} and q ^ {k + t + I, . . . ,mn}. 

On the other hand, since Wk{(f){{Eii + Ejj) ® In)) = [0, 2n/k], there exists a unitary matrix V of 
the form V = Is®Vi® Ik-t such that Yl^=i hp{V,j) = n, which implies 

hp{j) > hq{j) for all p G {1, . . . , s} and g G {/c + t + 1, . . . , mn}. 

Hence, we have 

^i(j) = ■ ■ ■ = hs{j) = /ifc+t+i(i) = • • • = /imn(i) for j = 7 + 1, . . . , m. 

Interchanging the roles of {2, . . . , 7} and {7 + 1, . . . , 27 — 1} and applying the same argument, we 
have ([I5|). 

Let r = + 1,. .. ,mn - 1}. Note that hi{U, 1) = ai(l) > a.^ni'^) = hmn{U, 1). We have 

mn 

J2hp{U,l) > J2 hp{U,l) = 0. 

peT p=k+l 

By the fact that hi{U,j) = hi{j) = /im„(j) = KmiUJ) for j = 2, . . . ,m, we have 

(m \ ml mn 

hp{u,i) + Y,hp{u,j)\ = Y,hp{u,i)+Y,{ 
i=7 / peT j=7 yp=fc+i 

= hp{U, l) + k>k 
which contradicts with Wk{(t>{{Eii + -£7,7 + • • • + Emm) ® In)) = 1- Hence, we get (fT2|) . 



PRESERVING THE HIGHER NUMERICAL RANGES 



13 



Next, suppose there is some 2 < j < m such that 4>{Ejj 
assume 4>{Eu <8> In) = Diag (ai(l), • • • ,amn.(l)) with 



In) ^ 0, say, j = 2. Again, we can 



ai(l) = • • • = afc+t(l) > afc+t+i(l) > • • • > amn(l), I <t <k-l 

and (j){E22 <^ In) = Diag (ai(2), . . . , amn(2)) with 

ai(2) > • • • > a,(2) > a,+i(2) = • • • = a^„(2) = 0, 1 < s < A;. 

Recall that there is a unitary matrix U satisfying (a), (b), and (c). Suppose Ui,U2 are unitary 
matrices such that U1B21UI = Diag (ai(2), . . . , afc(2)) and U2B12U2 = Diag (afc+i(l), . . . , am„(l)). 
Replace f/ with [Ui © U2)U. Then the new matrix U also satisfies (a), (b), and (c). Moreover, 



U(l)iEn^In)U* 

U4>{E22®In)U* 



Diag(ai(l),... 

(1)), 

Diag (ai(2), . . . ,Omn(2)) 



which implies that U is of the form U = (B (B with U3 € Mg, U5 E Mk_t. We can assume 
[/s = and C/5 = /fc_t. 

For any given j € {7 + 1, . . . ,m}, we denote by ai > • • • > the eigenvalues of Bji and 
/3i > • • • > the eigenvalues of Bj2. Then /3fc > ai- By Wk{(t>{{E22 + Ejj)®In)) = [0, 2n/k], there 
is a unitary V such that 



(18) 



y(S2i © B22)V* = y © and V{Bji © Bj2)V* = ^1 © ^2 



with y, Zi € Mfc, tr (y) = tr {Zi) = n. Suppose is a unitary matrix such that WYW* = B21. 
Replace V with (W © Ik)y ■ Then we still have (jlSp with Y = B21. Moreover, V is of the form 
= Vi © V2 with Vi S Ms and we can assume Vi = /g. Partition Bji and Zi as 



B 



'ji 



Cii C12 
C21 C22 



Zi 



Du D12 
D21 D22 



with Cii, Dii G Ms- We can rewrite the second equation in (jlSh as 



Cii (^12 
C21 (^22 



B 



D21 D22 



Z2 



It is clear that Du = Cu. Since tr (Cn + D22) = n equals to the sum of the k largest eigenvalues 
of Bji (BBj2, we see that tr (-D22) = Sp=i /5p- Applying Lemma [2T| we have D12 = -D21 — 0> which 



a. 



/3, 



k-s+l 



implies C12 = C21 = and a{Cii) = {ai . . . , as}- It follows that ai = ■ ■ 
/3k and Cn = ai/g. 

Similarly, considering IVfc(</'((^ii + Ejj) (g) /„)) = [0, 2n/k], there is a unitary V such that 

(19) y(Sii © Bi2)V* = ai(l)4 © y and © S^^)!^* = ^1 © ^2 

with Zi G Mfc, tr {Zi) = n. Suppose W is a, unitary matrix such that WYW* = B12. Replace V 
with (4 © W)V. Then we still have ([l9]) with Y = B^. Moreover, V is of the form V = ViBV2 
with V2 G Mk-t and we can assume V2 = Ik-t- Partition Bj2 and Z2 as 



B 



'J2 



R2I R22 



Zo 



Sll Si2 

S21 S22 
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with Rii, Sii € Mt. We can rewrite the second equation in (fT9]) as 
"Vi 

h-t 

Since tr (Zi) = n is the sum of the k largest eigenvalues of Bji © Bj2, which equals to the sum of 
the k largest eigenvalues of Cn © C22 © Rii, we see that the eigenvalues of Rn are also the t largest 
eigenvalues of Bj2- Hence, we have R12 = R21 = and i?22 = hik-t = ciilk-t- 
So we have 

(20) hi {]) = hi {U, j) = hmn{U, j) = hmn{j) for j = 7 + 1, . . . , m. 

Similarly, (pOj) holds for j = 2, . . . , 7. 

Again, we have pT|) . which contradicts with Wk{(l){{Eii+E^^^ + - ■ ■ + Emm)®In)) = 1- Therefore, 
(p{Ejj (8> In) has a negative eigenvalue for every I < j < m. 

Claim 5 . For any 1 < i < m and 1 <: j < n, the largest eigenvalue of (l){Eii (X" Ejj) is 1/k and, 
hence, ^Imn - (piEu ® Ejj) > 0. 

Given any 1 < i < m, by the previous claims, we can assume 

(l){Eii (g) /„) = Diag (oi (i), . . . , amn{i)) 

with ai{i) = ■■■ = ttk+iii) > ■■ > amn{i) and arnn{i) < 0. Denote by di{i,j), . . . ,dmnii,j) the 
diagonal entries of (j){Eii (8) Ejj). Then 

du{i, 1) = • • • = du{i, n) = 1 

for any T C {1, . . . ,k + 1} with \T\ = k. It follows that du{i,j) = 1/k for all 1 < u < k + 1 and 
1 < J < f^- 

Applying Lemma |2.H each <j){Eii (8) -Ejj) is of the form 

Diag {di{i,j),... ,4+i(«,i)) ® X{iJ), 
where the largest eigenvalue of X{i,j) is less than or equal to 1/k. Thus, we get the claim. 

Now, let ipiA (S)B) = (tr {A (g) B)/k)Imn - <P{A. B). Then Wk{ij{A B)) = Wk{A <S) B) for 
all A G Hm, B G Hn, and ilj{Eii iS) Ejj) > for all 1 < i < m, 1 < j < n. Applying the same 
arguments as in the first case on ip, we conclude that satisfies (1) and, hence, (j) satisfies (2). 
The proof is completed. □ 



Cu 



G 



22 



Vi 



Zi 



Rii R12 
R21 R22 



5*11 512 

5*21 »S'22 
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